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Abstract 

We report a measurement of the branching fraction for t“ —> 7r“7r*^ z/,- and the invariant mass 
spectrum of the resulting vr^vr*^ system using 72.2 fb“^ of data recorded by the Belle detector at 
the KEKB e"*“e“ collider. The branching fraction obtained is (25.15 ± 0.04 ± 0.31)%, where the 
first error is statistical and the second is systematic. The unfolded Tr~7r^ mass spectrum is used to 
determine resonance parameters for the p{770), /9'(1450), and p"{17QU) mesons. We also use this 
spectrum to estimate the hadronic contribution to the anomalous magnetic moment of the muon. 

TAGS numbers: 13.40.Gp, 13.35.Dx, 14.60.Fg 
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INTRODUCTION 


Among the decay channels of the r lepton, t~ —> 7r“7r° Ur has the largest branching 
fraction. The decay is dominated by intermediate resonances and thus provides information 
on the properties of the p(770), p'(1450), and p"(1700) mesons and their interference. Since 
leptons do not participate in the strong interaction, hadronic r decays provide a clean 
environment for studying the dynamics of hadronic states in an interesting energy range 
dominated by resonances. 

Under the Conserved Vector Current (CVC) theorem, the 7r~7r^ mass spectrum in this 
range can be used to improve the theoretical error on the anomalous magnetic moment of 
the muon = {g^j, — 2)/2. A recent review of the calculations of is given in Ref. P|. 

It is known that the theoretical error on is dominated by the (leading-order) hadronic 
contribution , given by the hadronic vacuum polarization. This contribution cannot 

be evaluated within the framework of perturbative QCD; however, it can be evaluated from a 
measurement of the cross section for annihilation to hadrons Si- Alternatively, CVC 
relates the properties of the tt+tt” system produced in —>• tt+tt” to those of the 7r“7r° 

system produced in r“—>7r“7r° decay; thus, using CVC and correcting for isospin-violating 

effects, r data have also been used to obtain a more precise prediction for ahad,LO _ 
Recently, data on —s-tt+tt" has become available from the CMD-2, KLOE, and SND 
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experiments M 0,0 @| • Data on r decays is available from the ALEPH 0,0, CLEO 
and OPAL [l^ experiments. The most recent evaluation of the hadronic contribution to 
using e+e“ data gives 0 = (25.2 ± 9.2) x 10“^°, while that using the r 

lepton data where applicable gives = (9.4 ± 10.5) x 10“^°. The experimental 

value is dominated by the BNL E821 measurement [0 (11 659 208 ± 5.8) x 10“^°. 
These differences correspond to 2.7 and 0.9 standard deviations, respectively, and thus there 
is a signihcant difference between the e’''e“-based and r-based predictions. To clarify the 
situation, more data for e’''e“ —> and for t~ —>■ 7r~7r^i>r decays are needed. In this 

paper we present a high-statistics measurement of the 7r“7r° mass spectrum produced in 


■71 71 ^ 


decays [1^ using data collected by the Belle experiment at the KERB 


e' e 


collider operating at a center-of-mass (CM) energy of 10.6 GeV. The data sample is about 
50 times larger than those of previous experiments. 


BASIC FORMULAE 


The differential decay rate for r —> vr 7r° z/,- can be expressed as 

l2 r , ir< ■. / \ 2 


1 dP _ _ 0 \ 67r|14rf| Sew 


Be 


. rrii, 




V-[S] 


( 1 ) 


where s is the invariant-mass-squared of the 7r“7r° system, n-(s) is the vector spectral func¬ 
tion characterizing the 7r“7r° system, \Vud\ denotes the CKM mixing matrix element, and 
Sew accounts for electroweak radiative corrections. Be and Bt^-k are the branching fractions 
for T~ —>• e~i'rVe and r“— z/,-, respectively. 

The corresponding 7r+7r“spectral function no(s) can be obtained from the e’''e“ —>■ tt+tt" 
cross section: 

a{e^e~ —*■ 7r’''7r“) = no(g), (2) 
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where s is the e'''e“ CM energy squared and oq is the hne-structure constant at s = 0. Up 
to isospin-violating effects, CVC allows one to relate the spectral function from r decays to 
the isovector part of the e+e" spectral function [l^ : 


n_(s) 



( 3 ) 


The mass spectrum of the two-pion system is typically expressed in terms of pion form 
factors; these are useful for comparing resonance shapes in the charged and neutral two-pion 
systems. The spectral function Vj(s) (j = —, 0) is related to the form factor F^(s) via 

^j(^) = ( 4 ) 

where /3-(s) (/3o(5)) is the pion velocity in the 7r“7r° (7r'''7r“) rest-system. The velocities /3j(s) 
are explicitly given by /3^(s) = [1 — (m.^- — /s] [1 — (m^- -|-m^o)^/s] and /3q(s) = 

[1 - 4m^-/s] . 

The leading-order hadronic contribution to the muon anomalous magnetic moment 
(^had,LO) jg the e’''e“ annihilation cross section via the dispersion integral 


hadjLO 


/oo^Y 

[ Stt J 




K{s) ds, 


( 5 ) 


where s is the invariant-mass-squared of the two-pion system, and R{s) = a{e'^e~ —*• 
hadrons)/(47raQ/3s). The kernel is a smooth function increasing from 0.63 at the 

threshold s = 4m^ to unity at s = oo |l7|. Due to the 1/s^ dependence, hadronic final states 
at low energy dominate the contribution to in fact about 70% of jg 

the two-pion state having 4m^ < s < 0.8 (GeV/c^)^. Consequently, the 27r spectral function 
in r data is useful to obtain predictions for Using Eqs. © and © to evaluate © 

we obtain 


a 


TTTT 


V 37r y 



3f(s) 


K{s) ds , 


( 6 ) 


where ”...” indicates the integral above the region. To determine V-{s), one must measure 
both the branching fraction for Vr and the 7r“7r° mass spectrum {l/N){dN/ds). 

Here we report new measurements for both of these. 


DATA SAMPLE AND SELECTION CRITERIA 

The data sample used was collected by the Belle detector at the KEKB energy-asymmetric 
collider [l8|. It is based on an integrated luminosity of 72.2 fb“^ recorded at a CM en¬ 
ergy of 10.6 GeV. The Belle detector is a large-solid-angle magnetic spectrometer consisting 
of a three-layer silicon-vertex detector (SVD), a 50-layer central drift chamber (CDC) for 
charged particle tracking and specific ionization measurement {dE/dx), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-ffight scintillation 
counters (TOE), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to identify muons and to detect 
mesons (KLM). The detector is described in detail elsewhere |l9|. 
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To study backgrounds and determine selection criteria, we perform Monte Carlo (MC) 
simulation studies for various processes. The KORALB/TAUOLA program BIHI is used 
for r’'‘r“-pair generation, the QQ generator for BB and qq continuum processes, the 
BHLUMI program for radiative Bhabha events, the KKMC program for radiative 
//■•■/i^-pair events, and the AAFHB program for two-photon processes. The BHLUMI 
and KKMC programs include higher-order radiative corrections and are among the most 
accurate programs available. The detector response is simulated by a GEANT-based pro¬ 
gram [^. In order to realistically simulate beam-induced background, detector hits taken 
from randomly-triggered data are added to wire hits in the CDC and to energy deposits in 
the ECL. 


r'*'r pair selection 

The event selection consists of two steps. Initially, a sample of generic e+e“ —*■ 
events are selected with relatively loose criteria. From this sample r “—Ur decays are 
then selected. The number of generic r+r“ events is used to determine the t~—>■ 7r~ 
branching fraction. 

Generic events are selected by requiring that the number of charged tracks in an 

event be two or four with zero net charge; that each track have a momentum transverse 
to the beam axis (pt) of greater than 0.1 GeV/c; and that each track extrapolate to the 
interaction point (IP) within ±1 cm transversely and within ±5 cm along the beam direction. 
To suppress background from Bhabha and events, the reconstructed CM energies and 

the sum of the momenta of the two leading tracks are required to be less than 9.0 GeV/c. 
The maximum pr among the tracks is required to be greater than 0.5 GeV/c. Beam-related 
background is rejected by requiring that the position of the reconstructed event vertex be 
less than 0.5 cm from the IP in the transverse direction and less than 2.5 cm from the 
IP along the beam direction. The polar angle of the leading particle with respect to the 
beam axis {6*) in the CM frame is required to be in the hducial region of the detector: 
35° < r < 145°. 

To reduce remaining background from Bhabha, pi ^and two-photon events, a cut is 
applied in the two-dimensional plane of the missing-mass MM and the direction of missing 
momentum in CM where MM is evaluated from the four-momenta of the measured 
tracks and photons: (MM)^ = (Pini — Ptracks — P-ysY- In this expression Pini is the four- 
momentum of the initial e’''e“ system. Each photon (reconstructed from clusters in the 
calorimeter) must be separated from the nearest track projection by at least 20 cm and have 
an energy greater than 0.05 GeV in the central part (—0.63 < cos 6 * < 0.85), and 0.1 GeV in 
the endcap part (—0.90 < cos 6 ' < —0.62 and 0.85 < cos 6 * < 0.95). Photons measured at the 
detector edge are rejected. A scatterplot of MM vs. for data is shown in Fig.[T] In this 
plot, events at MM k, 0 are due to radiative Bhabha events and e’''e“ —> /i^/i“( 7 ), while 
events in the high-MM region are from two-photon processes. Events within the octagonal 
region are selected as r+r“candidates. 

Candidate events are divided into two hemispheres in the CM frame with respect to the 
highest momentum particle, and the remaining background from e’''e“ annihilation processes 
is suppressed by selecting events with low multiplicity as characterized by the quantity 
Xpart = (^^tr + ^ 7)1 X (^tr + ^^ 7 ) 2 , where ntr,j and are the numbers of tracks and photons 
in hemisphere j. We require Xpart < 25. Finally, in order to eliminate Bhabha events in 
which one or both electrons produce a shower in material near the interaction region, the 
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FIG. 1: Missing mass (MM) versus the polar angle for the direction of the missing momentum 
(^miss)- shows MC c+c" —> T^T~ eveuts, and the right plot shows the data. Events 

inside the octagonal region are selected as r'’'r~-pair candidates. 


acoplanarity angle ^ between the first and second highest momentum tracks is required to 
be ^ > 1 °, where ^ = || 0 i — 02 | — 7 r| is defined as the two-track acollinearity in azimuth. 

After applying all selection criteria, 22.71 x 10® r'''r“-pairs survive. The background is 
estimated using MC simulation. The dominant source is from continuum processes —> 

qq {q = u,d,s,c) and amounts to 5.5%. Backgrounds from Bhabha events, and 

two-photon n~) events are estimated to be 0.6%, 0.4%, and 0.8%, 

respectively. Other sources are found to be small. These background estimates are checked 
by comparing the number of events in control samples. The control samples for continuum, 
Bhabha-|- and two-photon processes are high multiplicity events having 25 < Xpart < 

30 or |MM| < 0.5 GeV/c^ or |MM| > 8.0 GeV/c^, respectively. The differences in event 
yields for these control samples and the MC predictions (5-10%) are included as systematic 
errors for the results discussed in latter sections. 

r“— >7r“7r® candidate selection 

Within the r’''r“-pair sample, r“—>• decays are reconstructed by requiring that 

there be both one charged track and one tt® in one hemisphere. The tt® candidate is selected 
based on the normalized invariant mass — m^o)/c’‘p,.y, where is the mass 

resolution of the 77 system. The value of ranges from 0.005 GeV to 0.008 GeV, depending 
on the TT® momentum and polar angle. Pairs of photons with |S'p,p.| < 9.0 are considered as 
TT® candidates. To keep beam-related background at a negligible level, we require that the 
CM momentum of the vr® be greater than 0.25 GeV/c and the photon CM energy be greater 
than 0.08 GeV. 

The distribution in the normalized di-photon invariant mass for the selected 
sample, where there are one charged track and one tt® candidate in one hemisphere, is 
shown in Fig. |21 The lower-side tail of the distribution is primarily due to rear and 
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(m„ - m/)/a^ 


FIG. 2; Normalized 77 invariant mass spectrum in the data(points) and the MC 

event (histogram), for the sample described in the text. The plotted data correspond to 6.1% of 
the full data used in this analysis. The arrows indicate the signal region — 6 < < 5 and the 

sideband regions 9 < < 7. The sideband regions are used to subtract fake- 7 r° background. 


transverse leakage of electromagnetic showers ont of the CsI(Tl) crystals and the conversion 
of the photons in the material located in front of the crystals. Good agreement between 
data and Monte-Carlo indicates that these effects are properly modeled by the Monte-Carlo 
simulation. We define the interval —6.0 < < 5.0 as the 7 r° signal region. Spurious 7 r° 

background is small and estimated from the sideband regions 7 < \Sryry\ < 9. To reduce feed- 
down background from multi- 7 r° decays such as t~ —> 7r~{nn^)ur ( n > 2 ), signal candidates 
(in a hemisphere) are rejected if there are additional 7 ’s in the same hemisphere with energy 
greater than 0.2 GeV. 

The 7 r“ 7 r° invariant-mass-squared (M^^o) spectrum is obtained assuming the pion mass 
for the charged track; it is shown in Fig. |21 along with the MC prediction. To improve the 
energy resolution of the 7 r°, a 7 r° mass constraint is imposed. The amount of a spurious 
7 r° background depends on the region, varying from 4% to 15%. (This is subtracted 

using sidebands.) The final sample contains 5.55 x 10® t~^ h~'K^ Ur candidates after 
the TT® background subtraction, where h~ denotes tt or K. This sample is 50 times larger 
than those of previous studies. 

Feed-down background arises mainly from multi-vr® modes such as t~ — >■ 'K~{n'K^)i'r 
(5.5%) and r —> iF“ 7 r®z/^ (0.48%). Including other modes, the total feed-down background 
is (6.0 ± 0.1)%. The error listed includes statistical uncertainty as well as the uncertainty 
in relevant branching fractions. Background from non-r processes is negligible, except that 
from continuum processes. The amount of continuum background is estimated from MC 



















FIG. 3; Invariant-mass-squared (M^^o) distribution for Vt- The solid circles with error 

bars represent the data, and the histogram represents MC simulation (signal-|-background). The 
open area shows the contribution from r“—> 7r“7r® 1 ^,-; the narrow cross-hatched area shows that 
from T~ ^the wide cross-hatched area shows that from t~ ^ h~ and the striped 

area shows that from qq continuum and other non-r processes. 


simulation to be (2.45 ± 0.05)%. The normalization of the continuum MC is checked using 
data in the high-mass region 

The MC simulation of r decays is based on the TAUOLA program j^. A small difference 
observed between data and MC in Fig.3 for > 2.0 (GeV/c^)^ is attributed to the 

p"(1700) resonance, which is not included in the current TAUOLA program. 


MEASUREMENT OF THE BRANCHING FRACTION 


Formula 


The branching fraction for r -^h 7r° Ur (Bh-Ko) is determined by dividing the signal yield 
A^/ittO by the number of selected generic r'''r“-pairs Nrr- 


Bhir^ 




(1-h 


feed-down 




(1 - brr) 



(7) 


In this formula, b^r is the background fraction in the sample, e^r is the efficiency of 

the r'''r“-pair selection, ej^^o is the efficiency for t~ —> h~'K^v decays to pass the r'*’r“-pair 
selection, and is the efficiency for t~ —*■ h“7r°z/ decays satisfying the r+r'-pair selection 
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to pass the h“7r° selection. The product is the overall detection efficiency for the 

h~TT^v final state. The parameter ft^eed-down jg fraction of h~7r^i> candidates coming from 
other r decay modes, and is the fraction coming from non-r processes. In this formula, 
several common uncertainties such as that in the luminosity, that in the cross section for 
T^T~-Y>aAi production, that in the trigger efficiency, and that in the selection efficiency 
cancel in the ratio. The value for each factor is listed in Table H] along with the statistical 
error. 


Parameter 


Values 


Ett 


h-K^ 

fb = 


Err 


30.81 ± 0.05 % 
34.26 ± 0.07 % 
1.112 ±0.003 


e 

b 

b 

b 


ID 


TT 

feed-down 

non-T 

/iTT^ 


42.62 ±0.13 % 
7.66 ± 0.03 % 
5.98 ± 0.08 % 
2.45 ± 0.06 % 


TABLE I: Values of parameters used for the branching fraction measurement along with statistical 
errors. 


Systematic uncertainty 

The sources of systematic uncertainty are listed in Table E] The uncertainty in the 
tracking efficiency is estimated using D*^ K~7r^7r^ decays to be 1% per track. 

A large part of this uncertainty cancels in the ratio of Eq. ©; the resulting uncertainty 
from this source is = 0.12 %. The detection efficiency is obtained from the 

ratio of > K~7r~^n^ to —>■ K~7r~^ decays, in which the branching fractions are 

known precisely. The uncertainty is estimated to be ±1.7% for a 7r° momentum less than 
1.0 GeV/c. As a consistency check, the branching fraction is re-measured after changing 
the photon threshold from 0.05 GeV to 0.10 GeV; the difference in is only 0.10%. 
The uncertainty in background in the non-r sample is estimated from the control 

sample above the r mass region, while the uncertainty in feed-down background 
is obtained from the uncertainty in r“ —h“(n7r°)z/^ and t~ ^ branching fractions. 

The veto of additional y’s is required in the event selection to reduce background from 
multi-7r° decay channels. However, it also vetoes signal if there are photons radiated in the 
initial or final state. In addition, photon candidates can also appear due to electromagnetic 
shower fragments and/or mis-reconstructed of electrons. The uncertainty from these sources 
is estimated by changing the veto threshold by ±0.1 GeV; the resulting change in is 
only ±0.05%. Signal events are flagged by several trigger conditions that require two or more 
GDG tracks with associated TOF hits, EGL clusters, or a significant sum of energy in the 
EGL. This redundancy allows one to monitor the efficiency of each trigger requirement. The 
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Source of uncertainty 

(%) 

Tracking efficiency 

0.12 

vr^/y efficiency 

0.25 

Background for 

0.09 

Eeed-down background for r“— 

0.04 

Non-r background for r“— u-y 

0.05 

7 veto 

0.05 

Trigger 

0.08 

MC statistics 

0.04 

Total 

0.31 


TABLE II: Systematic uncertainties for the r —> h -k^ v^- branching fraction. 

uncertainty arising from the trigger is estimated by assuming there is a ±3% uncertainty in 
the track and energy trigger efficiencies, which is the maximum variation measured during 
experimental running. The resulting uncertainty on is small (0.08%) since the 
trigger efficiency is high (97%). 


Results 

Inserting all values into Eq. dZD gives 

Bh^o = (25.60 ± 0.04 ± 0.31)%, (8) 

where the first error is statistical and the second is systematic. Subtracting the small 
kaon-channel branching fraction listed in the PDG j2^ [Bk--kO = (0.45 ± 0.03)%] gives a 
branching fraction of 

B^^o = (25.15 ± 0.04 ± 0.31)%. (9) 

This result is in good agreement with previous measurements, as shown in Table IIIII Our 
statistical error is significantly lower than those of the other measurements, while our sys¬ 
tematic error is similar to those of the others (except for the ALEPH result). 


Experiment 

BhA%) 

Reference 

OPAL 

25.89 ±0.17 ±0.29 

[ 22 ] 

ALEPH 

25.924 ± 0.097 ± 0.085 

[ 22 ] 

L3 

25.05 ±0.35 ±0.50 

[30] 

CLEO 

25.87 ±0.12 ±0.42 

m 

This work 

25.60 ±0.04 ±0.31 



TABLE III: 


Branching fractions for r —measured by various experiments. 
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FIG. 4: Invariant-mass-squared (M^^o) distribution for r —>7r I'j- after background subtrac¬ 

tion. 


MEASUREMENT OF THE MASS SPECTRUM 

In order to obtain the true 7r~7r^ mass spectrum, one must correct for (1) background, (2) 
smearing due to finite resolution and radiative effects, and (3) mass-dependent acceptance. 


Background Correction 

As noted earlier, background entering the r“—> 7r“7r° z/,- sample is small. The sidebands 
of the M-y-y distribution are used to estimate the fake 7r° contribution. This background 
dominates at values of less than about 0.25 (GeV/c^)^. 

As seen in Fig. El feed-down background arises from t~ —>■ h~{n7i^)h>T- and t~ — 
decays; both backgrounds dominate at low values of In the high mass region, con¬ 

tinuum background dominates. For this analysis we did not use information from particle 
identihcation (FID) detectors to separate charged pions from kaons, as the feed-down back¬ 
ground is dominated by r“ —> h~{n'n'^)ur rather than t~ —>■ K~n^Ur. The distribution 
after subtracting this background is shown in Fig. E) 


Acceptance Corrections 

The detector effects include M^^o-dependent acceptance and bin-by-bin migration caused 
by the hnite mass resolution. We correct for these effects by performing an unfolding proce- 
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FIG. 5: Correlation between Mgg^ (vertical axis) and (horizontal axis), the generated and 

observed invariant masses squared of the system in t“—> 7r“7r® decay. 


dure. The unfolding program used is that employed by the ALEPH experiment • In this 
program, the unfolding is based on the Singular-Value-Decomposition (SVD) method js^, 
in which the acceptance matrix is inverted by limiting the number of singular values to 
only those that are statistically signihcant. The output of the program is the unfolded 
distribution and its covariance matrix. 

The correlation between the generated quantity and the measured one is shown 
in Fig. El The hgure shows a clear correlation between and The resolution in 

is 0.005 (GeV/c^)^ in the low-mass region and 0.030 (GeV/c^)^ in the high-mass region; 
thus the bin size chosen is AM^ = 0.050 (GeV/c^)^ so that the off-diagonal components of 
the acceptance matrix are small. 

The acceptance as a function of is shown in Fig. El The acceptance varies smoothly 
and its average value is 17%. It decreases at low values of due to the overlap of 7 
clusters with the vr” track projection at the calorimeter. 


Results 

The unfolded s = ^ spectrum dN/ds is shown in Fig. [T] The square roots of the 

diagonal components of the covariance matrix are used as the errors. The p peak and the 
shoulder due to the p'(1450) are clearly visible. The dip at s ~ 2.5 (GeV/c^)^ is caused by 
destructive interference between the p'(1450) and p"(1700) resonances. 

To obtain parameters for the p, p' and p resonances, a Hi using Breit-Wigner functions 
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FIG. 6: Acceptance as a function of as determined from r —>7r tt^ I'j- MC simulation. 


is performed. Since the unfolded mass spectrum has bin-by-bin correlations, the off-diagonal 
components of the covariance matrix X are included in the evaluation; 

= Y1 iVi - «)) iVj - f{sj; a)) , ( 10 ) 


where yt is the measured value at the z-th bin, f{s]a) is the value of the function for 
parameters a, and is the inverse of the covariance matrix. 

In the £t, the s dependence of the decay rate is given by Eq. (HD- The pion form factor in 
Eq. is parametrized with Breit-Wigner functions corresponding to the p, p'(1450), and 
p"(1700) resonances: 

(BWp + (3.BW,,+^- BW ,.), ( 11 ) 

1 -I- p -1- 7 


where the parameters (3 and 7 (denoting the relative size of the two resonances) are in 
general complex. We use the Gounaris-Sakurai (GS) model for the Breit-Wigner shape: 


BW^^ 


+ d-M,r,(g) 

{Mf - s) + /(s) - iysri(s) ’ 


with an energy-dependent width 


r.(s) 




( 12 ) 


(13) 
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FIG. 7; Fully-corrected distribution for T~^Tr~7r^ Ut- The solid curve is the result of a fit 

to the Gounaris-Sakurai model with /?(770), /3'(1450), and p”{1700) resonances. All parameters are 
floated. 

Here, k{s) = is the pion momentum in the 7i~7i^ rest frame. The functions /(s) 

and h{s) are dehned as 


/(s) 

h{s) 


Mf 

k%Ml) 


k‘^{s) {h{s) 


2 kis) y/s + 2k(s) 

- -IT) - 

71 y/s 2m^ 




ds 


5=M2 


(14) 

(15) 


with dh/ds\j ^2 = h{Mf) ( 8 fc^(Mf)) 


-1 




+ (27rMf) ^ and 


^ 3 ml M, + 2k{Mf) M, 

7T k^{Mf) 2m^ 27Tk{Mf) 7ik'^{M‘l) 

There are ten parameters in this formula: the masses (Mj) and the widths (Tj) for the p, p', 
and p” resonances, their relative amplitudes \(3\, I 7 I, and their phases (j)^ and 

All parameters are floated in the £t. When evaluating the the 1% systematic un¬ 
certainty resulting from the unfolding procedure is included in the diagonal part of the 
covariance matrix. This uncertainty is estimated by applying the same unfolding procedure 
to MC events and comparing the unfolded spectrum with the original. The result of the 
£t is shown in Fig. [7|as the solid line; the values obtained for the parameters are listed in 
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Parameter 

Fit result 
(all free) 

Fit result 
(fixed py) 

Mp (MeV/c2) 

774.6 ±0.2 ±0.3 

774.3 ±0.2 ±0.3 

Tp (MeV) 

150.6 ±0.3 ±0.5 

150.0 ±0.3 ±0.5 

Mp, (MeV/c2) 

1336 ± 12 ± 23 

1436 ± 15 ± 23 

Pp (MeV) 

471 ± 29 ± 21 

553 ± 31 ± 21 

m 

0.090 ± 0.009 ± 0.013 

0.161 ±0.020 ±0.013 

pp (degree) 

123.7 ±5.0 ±7.0 

149.1 ±2.4 ±7.0 

Mpn (MeV/c2) 

1600 ± 13 ± 4 

1804 ± 16 ± 4 

Tpn (MeV) 

255 ± 19 ± 79 

567 ± 81 ± 79 

ItI 

0.062 ±0.015 ±0.015 

0.136 ±0.024 ±0.015 

(degree) 

-64.1 ± 7.9 

[0] 

xVd-o.f 

55/51 

94/52 


TABLE IV: Results of fitting the distribution for r“ —> i/,- to the Gounaris-Sakurai 

model with the p{770), /9'(1450), and /9"(1700) resonances. The results for two cases, all parameters 
floated (the second column) and fixed cj)^ (the third column) are shown. For both cases, the first 
error is statistical and the second one is systematic. The systematic errors include the uncertainty 
of the backgrounds, unfolding, as well as the uncertainty of the photon energy scale. 


Table IIVI The results are compared with the previous ALEPH measurements in Table 0 
In the table, the hrst error is statistical and the second one is systematic. 

The main sources of systematic uncertainty are the photon energy scale, the unfolding 
procedure, and the background subtraction. The uncertainty in the p mass (0.3 MeV) is 
mainly due to the uncertainty in the photon energy scale. The uncertainty in background 
dominates for the p" parameters. Our result for the mass of the p resonance agrees well with 
the ALEPH and CLEO^^ results. 

As can be seen from the second and third columns of TablelVl where the interference angle 
(f)^ or three parameters Mpnpn are fixed, respectively, as in the previous ALEPH [23 
fit, the values for p' and p" resonance parameters are quite sensitive to the values of other 
parameters fixed in the fit. 

The results are shown in terms of the pion form factor squared (|iA(s)p) in Figs. |H1 
andini A dip caused by destructive interference between the p'(1450) and p"(1700) is clearly 
visible. For the first time production of the p'''(1700) in t~ decays has been unambiguously 
demonstrated and its parameters determined. For comparison, the figures also show results 
from the CLEO and ALEPH experiments, respectively; there is good agreement 
with both data sets. Figure ^1 shows our data and that of CLEO for the mass range 
0.2-2.2 GeV/c^, where the contribution to jg laj-gegp 
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Parameter 

This work 
(fixed 4>j) 

This work 
(fixed Mp/,, Tp „, 

ALEPH(t) 

Mp- (MeV/c^) 

774.3 ±0.2 

773.9 ±0.1 

775.5 ±0.7 

Pp- (MeV) 

150.0 ±0.3 

150.8 ±0.3 

149.0 ±1.2 

Mp, (MeV/c2) 

1436 ±15 

1395 ±4 

1328 ±15 

Pp' (MeV) 

553 ±31 

411 ±9 

468 ±41 

l/3| 

0.161 ±0.020 

0.095 ±0.02 

0.120 ±0.008 

(t)fj (degree) 

149.1 ±2.4 

161 ±2.0 

153 ±7 

Mpn (MeV/c2) 

1804 ± 16 

[1713] 

[1713] 

Pp// (MeV) 

567 ±81 

[235] 

[235] 

I 7 I 

0.136 ±0.024 

0.045 ± 0.002 

0.023 ±0.008 

(j)-^ (degree) 

[0] 

[0] 

[0] 

^^/(d.o.f) 

Reference 

94 / 52 

134/54 

119 / 110 
[29] 


TABLE V: Comparison of our fit results for the /9(770), /o'(1450), and /o'^(1700) parameters with 
those obtained by the ALEPH experiment [ 2 ^. The numbers in brackets indicate the values fixed 
in the fit. 


EVALUATION OF 


Using the unfolded s = ) distribution {l/N){dN/ds), the spectral function v is 

obtained by taking the inverse of Eq. (fT|l : 




ml 

( B,r7r\ 

6vr WrfPS'Ew 

1 Be ) 




2s \ 1 dN 

) N ds 


(U) 


The resulting function can be inserted into Eq. & to obtain the dominant low-mass con¬ 
tribution to the hadronic part of the anomalous magnetic moment, This assumes the 
CVC relation Q holds. 

There are several external parameters in these equations; the values used for these are 
listed in Table EH For nir, Vud, and Be, PDG values are used. For the electroweak 
radiative correction Sew, we use the recent value 1.0233 ± 0.0006, which is based on a 
consistent treatment of the isospin-breaking correction j3, 13. For the 7i~7i^ branching 
fraction, our measurement is consistent with the world average given in Ref. na. 

Including our result and the recent ALEPH measurement, the new world average is 


B^^o = (25.42 ±0.11)%. 


( 18 ) 


We use this new world average for the evaluation of 

The integration in Eq. is carried out numerically by taking the sum of the integrand 
evaluated at the center of each bin. The statistical error in is calculated including the 
off-diagonal elements of the covariance matrix Xij: 





( 19 ) 
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FIG. 8: Pion form factor for t“—> 7r“7r^ Vr- The solid circles show the Belle result and the open 
squares show the CLEO result The error bars for the Belle data include both statistical and 
systematic errors added in quadrature. The solid curve is the result of a fit to the Gounaris-Sakurai 
model with the /?(770), /9'(1450), and p"(1700) resonances, where all parameters are floated. 


Source 

Value 

Relative error 

(%) 

AaY 

(10-10) 

Reference 

Sew 

1.0233 ± 0.0006 

0.06 

± 0.32 

mm 

Vud 

0.9734 ± 0.0008 

0.08 

± 0.42 

mn 

Be 

(17.84 ± 0.06)% 

0.34 

± 1.82 

mn 


(25.42 ± 0.11)% 

0.43 

± 2.30 


Total external 



± 3.0 



TABLE VI: Values of the external parameters and the systematic errors for arising from these 
sources. 


Because of uncertainties associated with the background estimate and with the acceptance 
correction in the lowest mass region, the integration is carried out over the mass range a/s = 
0.50-1.80 GeV/c2. 
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FIG. 9: Pion form factor for t ^ tt tt^ Vr- The solid circles show the Belle result and the open 
squares show the ALEPH result 0 . The error bars for the Belle data include both statistical and 
systematic errors added in quadrature. The solid curve is the result of a fit to the Gounaris-Sakurai 
model, where all parameters are floated. 


Systematic uncertainty 

Systematic uncertainty in arises from both external and internal sources. The errors 
arising from external parameters are summarized in Table IVll the total systematic error 
from these sources is ±3.0 x 10“^° (dominated by 

The systematic error arising from internal sources (specihc to this measurement) are 
listed in Table IVTTl There are two sources of background in the vr 7r° sample: (i) feed-down 
from T~ —> h~{n7r^)i>r and t~ — K~7r^i>r, and (ii) non-r background. In the hrst case, 
the uncertainty in the branching fraction is used to estimate the error. In the second case, 
the uncertainty in the background as estimated from the control samples is assigned as the 
error. As mentioned earlier, the fake-7r° background is subtracted using sideband events; 
the uncertainty is determined by varying the signal and sideband regions. 

The ratio of the branching fractions for the decays —>■ and > K~ti^ 

is used to monitor the vr® efficiency. It is found that the shape of the mass spectrum is 
insensitive to uncertainty in the vr® efficiency, as it is only at the few % level. Adding all 
individual errors in quadrature gives a total error on arising from internal sources of 
±1.0 X 10“^®. 

To check the stability of we perform the following tests: 

1. The sample is divided into subsamples based on the tag-side topology, i.e., one electron, 
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FIG. 10: Comparison of the pmn form factor squared |Fn-(s)p measured by Belle to that measured 
by CLEG Q and ALEPH 
from the fit of the Belle t~- 


29| experiments in the p(770) and /9'(1430) mass region. Difference 


•vr TT*^ 


Vr data divided by the fit value is plotted. 


one-prong, or three-prong. The values of obtained from these subsamples are 
consistent within the statistical errors. 

2. The sample is divided into subsamples based on the running period, e.g., years 2000, 
2001, or 2002. Again, the values of obtained are consistent within the statistical 
errors. 

3. The sample might be sensitive to the requirement on the overlap region between the 
projection of the charged track and 7 clusters. To estimate this sensitivity, we select 
events with a tighter isolation requirement on 7 ’s and on the track extrapolation; 
50 cm instead of 20 cm. 

The resulting variation in is small and is included as an additional systematic 
error. 


RESULTS 

The result for integrated over the mass range ^/s =0.50-1.80 GeV/c^ is 

a;j^[0.50,1.80] = (464.4 ± 0.6 (stat.) ± 1.0 (sys.) ± 3.0 (sys. ext.) x 10"^°, 

where the hrst error is statistical and the second and third errors are systematic errors 
arising from internal and external sources, respectively. In addition, there is a systematic 
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Source 

X 10^° 

(0.50-1.80 GeV/c^) 

Background: 


non-T (e'''e“ — > qq) 

±0.11 

feed-down /i(n7r*^)n 

±0.09 

feed-down K~Tr^u 

±0.15 

Energy scale 

±0.10 

TT^/y selection 

±0.24 

7 veto 

±0.93 

Efficiency: 


7r°/7 

±0.35 

charged track 

< 0.10 

Integration procedure 

< 0.10 

Total internal 

± 1.04 


TABLE VII: Systematic errors for arising from internal sources (specific to this measurement). 


uncertainty caused by isospin violation effects arising from p-uj interference, from the vr 
and 7r° mass difference, and from radiative corrections (see Ref. j3). The overall correction 


is estimated to be (—1.8 ± 2.3) x 10“^°, where the central value is taken from Ref. 11| and 
we enlarged the error according to the value in Table 5 of Ref. jH ; this correction is small 
because the threshold region is not included. Applying this correction gives 


a7[0.50,1.80] = (462.6 ± 0.6 (stat.) ± 3.2 (sys.) ± 2.3 (isospin)) x 10"^° , 

where the hrst error is statistical, the second is systematic, and the third arises from isospin 
violation. 

This result can be compared to those from previous r and experiments js^: 

a^’^[0.50,1.80] = (464.0 ± 3.0 (exp.) ± 2.3 (isospin)) x 10“^° (r : ALEPH, CLEO) 
a^^[0.50,1.80] = (448.3 ± 4.1 (exp.) ± 1.6 (rad.)) x 10"^° (e+e” : CMD2, KLOE). 

The hrst error includes both statistical and experimental systematic errors added in quadra¬ 
ture. The second error in the e’''e“ result is due to radiative corrections. Our result agrees 
well with the r-based result but is noticeably higher than the e’''e“ result. This supports the 
hypothesis that there is a difference between the mass spectra of the 27r systems produced 
in r-decay and —>■ tt+tt” reactions. 


In summary, we have studied the decay r“—>• 7r“7r°z/,- using high statistics data taken 
with the Belle detector at the KERB collider. The branching fraction is measured with 
1.2% accuracy, which is better than that in the previous experiments (except for the ALEPH 
result). In the unfolded 7r“7r° mass spectrum, in addition to the /o(770) and p'(1450) mesons. 
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the production of the p"(1700) in r“ decays has been unambiguously demonstrated and its 
parameters determined. The unfolded spectrum is used to evaluate the 2tt contribution to 
the muon anomalous magnetic moment in the region ^/s = 0.50 — 1.80 GeV/c^. Our 
results agree well with the previous r based results but are higher than the e+e“ results. 
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